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Environment induced decoherence entails the absence of quantum interference phenomena from 
the macroworld. The loss of coherence between superposed wave packets is a dynamical process the 
speed of which depends on the packet separation: The farther the packets are apart, the faster they 
decohere. The precise temporal course depends on the relative size of the time scales of decoherence 
and other processes taking place in the open system and its environment. We use the exactly 
solvable model of an harmonic oscillator coupled to a bath of harmonic oscillators to illustrate 
various decoherence scenarios: These range from exponential golden-rule decay for microscopic 
superpositions, system-specific decay for larger separations in a crossover regime, and finally the 
universal interaction-dominated decoherence for ever more macroscopic superpositions investigated 
in great generality in the accompanying paper [W. T. Strunz, F. Haake, and D. Braun; Phys. Rev. 
A, henceforth referred to as [SHB]]. 
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I. INTRODUCTION 

A collection of N + 1 harmonic oscillators N of which 
are mutually free but all coupled, symmetrically and har- 
monically, to the remaining oscillator, enjoys consider- 
able popularity as a model of an A^-freedom environment 
acting on a single- freedom system HJ, |^]. The impor- 
tance of the model lies in its rigorous explicit tractabil- 
ity through a normal-mode analysis. Most applications 
aim at revealing the effectively irreversible behavior of 
the central oscillator brought about by the coupling to 
the environment (alias "bath") for large N; in the limit 
N — > oo even strict irreversibility results under certain 
assumptions for the distribution of bath frequencies, and 
the central oscillator becomes linearly damped. 

While obliged to the tradition just pointed to, the 
present paper picks up a more recent trend and exploits 
the rigorous tractability of the model to reveal the emer- 
gence of classical behavior from quantum dynamics in the 
passage from the microscopic to the macroscopic world 
|^|, |], In particular, we study the temporal fate of su- 
perpositions of two wave packets for the central oscillator 
as in jj| , yet concentrate on the various decoherence sce- 
narios emerging as the distance of the superposed wave 
packets is varied from microscopic to macroscopic scales. 
As is by now well known such superposed packets loose 
their relative coherence, due to the dissipative influence 
of the bath, the faster the larger the initial separation of 
the two wave packets. The life time Td ec of the relative co- 
herence is inversely proportional to a power of the initial 
separation, Td cc oc (\/d) v with v > 0; within the power 
law, the separation d is referred to a microscopic quan- 
tum scale A, and therefore exceedingly rapid decoherence 
results as d is increased towards macroscopic magnitude. 
Keeping all other parameters of the problem fixed, it is 
the initial separation between the two wave packets only 
that determines the size of the decoherence time scale 
Tdcc relative to other relevant system or environmental 



time scales. 

The most familiar golden-rule limit 

7"sys 7"dec < ^~~ Tdiss ) 



(1.1) 



which allows the system to undergo many cycles during 
decoherence, can hold only as long as the separation be- 
tween the two wave packets remains below a certain limit. 
On increasing the separation we encounter a qualitatively 
different regime in which decoherence is faster than any 
system time scale, 



7"dec ^~sys 5 7"diss 5 



(1.2) 
In that 



irrespective of the relative size of r sys and r^h 
interaction dominated limit the free-motion Hamiltonian 
of the central oscillator, rather than the interaction with 
the bath, behaves like a weak perturbation during deco- 
herence. For yet larger separations decoherence becomes 
the fastest process by far, faster even than reservoir time 
scales, 



7"dec < ^ Tres ? 7"sys ) 7"diss • 



(1.3) 



The universal behavior resulting in the limit (1.3) may 
look like instantaneous decoherence on the classical time 



scales r B - 



, 1~di! 



While certainly requir ing separations d 



huge on the quantum scale A, the limit (1.3) will turn out 
to allow, surprisingly, moderate or even small d relative 
to every-day macroscopic scales. 

The pure limiting cases mentioned above allow for an- 
alytical treatment for general open systems, as shown in 
the accompanying paper [SHB]. Crucially, decoherence 
in the interaction dominated cases (1.2) or ( |1.3| ) becomes 
independent of the system Hamiltonian and may thus be 
regarded as the origin of the universally observed absence 
of quantum interferences in the macroworld. All of these 
general findings of the accompanying paper will be illus- 
trated for the oscillator model here. However, the prin- 
cipal purpose of the present paper is to study the inter- 
esting crossovers between the three regimes mentioned; 
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since these crossover effects elude the general asymptotic 
methods of [BHS], the exact tractability of the oscilla- 
tor model allows precious insights into the emergence of 
classical behavior. 



II. DECOHERENCE OF SUPERPOSED WAVE 
PACKETS 

We illustrate the various decoherence scenarios for a 
superposition of Gaussian wave packets 



\<P) = Ci\<Pi) +C2IV2) , |ci 
In the position representation, 
1 



C2 



1. 



(q\ipi) = tpi(q) 



(2 7 ra) 1 /4 



iPi(l-9i)/?> p -(g-9i) /4<r 



(2.1) 



(2.2) 



with i = 1,2. These packets are located in position 
space at qi with (rms) uncertainty Aq = ^fa and in 
momentum space at pi with uncertainty Ap = h/2^/a. 
We choose coherent states |j) with the minimum un- 
certainty AgAp = h/2 such that both Aq and Ap are 
Vn. To ensure good separation we stipulate that ei- 
ther Aq <C \qi — (fel or Ap <C |pi — p 2 | or both (see Fig. 
1). 



iq,- qj 





position 



FIG. 1: Position space density |i/>(g)| 2 for a superposition of 
two Gaussian wave packets with mutual distance much larger 
than individual spread. 



Th e initial density operator corresponding to the state 
(2.1) is a sum of four terms, 



Psys(O) 
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c i c* j \v i ){vA=Y, c i c *3$yM> ( 2 - 3 ) 



two "diagonal" ones weighted by probabilities \ci\ 2 
and two off-diagonal "interference terms" Psys(O) = 



IVi)(V2| 
and c*C2- 



Psy S (0) t weighted by the "coherences" c\c\ 



As in [SHB] we may employ the norm 



(2.4) 



as an indicator of the temporal fate of the relative coher- 
ence between the two superposed wave packets. Clearly, 
if the system were closed its unitary time evolution would 
leave that norm constant in time, N±2(t) = 1; interaction 
with a many-freedom environment will cause decay. 



III. HARMONIC-OSCILLATOR MODEL 

Our aim is to illustrate various decoherence scenarios 
for quantum superpositions \ip) — ci\ipi) -\-c2\f\) of wave 
packets of a harmonic oscillator of mass M and frequency 
Vt. Position and momentum operators Q and P obey the 
usual commutation relation [Q, P] — ih. The reservoir 
is a collection of harmonic oscillators as well, the ith of 
which has the coordinate Qi and the momentum P^, fre- 
quency u>i, and mass m; the coupling is taken bilinear in 
the positions, such that the three terms in the Hamilto- 
nian H — H sys + H rcs + H int read 



P 2 



1 



2M 2 



M0 2 Q 2 



*~ = E (& + \™4<X 



(3.1) 



flint = QB = Qj29iQi- 



This model and variants thereof have been used exten- 
sively over the years to investigate dissipative quantum 
dynamics jl], ^, ^ , and decoherence in particular || . Its 
popularity is due to the fact that it allows for an explicit 
exact solution for the many-body Schrodinger equation. 

The dissipative and decohering influence of the reser- 
voir is encoded in the thermal autocorrelation function 
of the bath coupling agent B, 



W)B(0)) = 



2muj, 



[(2nth(wj) + 1) cosujit 



1 sin uj. 



t], 



(3.2) 

where the time dependence in B(t) — e lHrost / h Be lH «= s */ n 
refers to the free motion of the bath and n t h(w) = 
(e hu] / kT — 1) _1 is the thermal number of quanta in an 
oscillator with angular frequency to. Assuming the num- 
ber N of bath oscillators to be large, it is customary j^] 
to introduce the spectral density 



/ ' * mi,): 



muji 



(3.3) 



such that the real and imaginary parts of the correlation 
may be expressed as 



{\{B{t),B}) = - 



du) J(uj)(2n t h(to) + 1) cosut 
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£[B(t),B]) = - 
n 7r 



cLuJ(lo) shout. 



(3.4) 



Note that the imaginary part (the damping kernel) is 
independent of h, while the real part, which describes 
equilibrium fluctuations, becomes independent of h only 
in the high temperature limit when kT/h is the largest 
frequency involved. For small temperatures, however, 
only quantum fluctuations remain, such that the real part 
of the correlation function becomes of first order in h. A 
so-called Ohmic bath is provided by a spectral density 
with a linear frequency dependence for small u>, 



J(lu) = MjloMlu/A) 



(3.5) 



with a cutoff function such that / c (0) = 1. The rate 7 is 
a measure of the coupling strength and turns out to be 
the classical damping rate; A is a cutoff frequency. For 
the model to be physically sensible, the cutoff frequency 
A is assumed much larger than the frequency f2 and the 
damping constant 7 such that 



f c (n/A) « 1 



(3.6) 



We shall use f c (x) = 1/(1 + x 2 ) 2 for our simu- 
lations. Our model gives an initial value (B 2 ) = 
^ j-~ dcuLu(2n th (Lu) + 1)/ C («/A). In the high- 
temperature limit, we obtain 

2 f°° 

(B 2 ) = MkTjA— / dxfJx) = MkTjA ■ 0(1) , (3.7) 

7T JO 

the remaining integral being a real number of order one. 
The zero-temperature limit 



(B 2 



= hMjA 2 - 

7T 



dxxf c (x) = HMjA 2 -0(1) (3. 



mainly differs from the high-temperature one by the re- 
placement of the thermal energy kT with the cutoff en- 
ergy fiA. The choice of frequencies 0,7, A, kT/h deter- 
mines the system and reservoir time scales of our model. 
We choose 7 = l(T 5 fi, A = 10 2 fl, kT/h = 20fi. 



A. Dynamics: Exact Master equation 

The dynamics of the system is determined by the re- 
duced density operator p sys (t) = pt- Assuming initial 
decorrelation of system and bath, we encounter the fol- 
lowing well known evolution equation fel 



Pt 



1 



[H sys ,p t ] 



(3.9) 



+ ||[^]H 

+%[Q,[P,Pt] 
h 



2i^ L 
_ dt 
ti 



-{P,Pt}] 



2 [QdQ,pt}} 



with real-valued time dependent functions |lC| ] 
a t ,b t ,c t ,d t whose physical meaning as drift coeffi- 
cients (a, b) and diffusion coefficients (c, d) will become 



clear presently; they approach constant values on the 
bath correlation time scale 1/A. 

Remarkably, despite the generally non-Markovian na- 
ture of the true open system dynamics, the evolution of 
the exact /? sys (i) is governed by a time-local differential 
equation; memory effects arc encoded in the time depen- 
dent coefficients at, ■ ■ ■ ,dt- 

We need not specify the precise time dependence of all 
four coefficients at this stage, but would like to mention 
their behavior at early times, 



h 

Ct 

d t 



0{t 2 ) 
0{t 3 ) 
1 

2M 



(B 2 )t- 



(3.10) 



2\ .2 



(B 2 )t 



0(t 4 



0(t 3 ) . 



It is this early time dependence of the (diffusion) coef- 
ficients Ct and dt that will turn out relevant for the de- 
coherence of the largest superpositions alias Schrodinger 
cat states, to be discussed later. 



B. Wigner Representation 

It is useful to switch to a phase-space representation 
and to express the above master equation (3.9) for p as 
an evolution equation for the Wigner function 



W(q,p) 



1 

2nh 



dv{q-u/2\p\q + y/2)e lup ' h . (3.11) 



For the init ial superposition of two Gaussian wave 
packets as in (2.1) the Wigner function has three dis- 
tinctive features as displayed in Fig. 2: two Gaussian 
wave packets in phase space arising from the diagonal 
terms p 11 and p 22 in ( |2.3D , and an oscillating pattern 
in between the two Gaussians, due to the coherences p 12 
and p 21 . Decoherence leads to the disappearance of those 
oscillations, measured nicely by the norm Ni2(t) as we 
will reveal shortly. 

In terms of the Wigner function, ( |3.9| ) takes the form 
of a Fokker-Planck equation 



Wt = 



d 2 d 2 \ 
-°t^-^r + d tTr^ w t . 



a t )q + b t p} 



(3.12) 



We now read off the meaning of the coefficients a, b, c, d 
in the evolution equations: while the term involving at is 
a mere potential renormalization due to the coupling, the 
term involving bt describes damping. The two remaining 
terms represent diffusion; there is a mixed second-order 
derivative with coefficient ct, while the time dependent 
momentum diffusion involving d t reflects the stochastic 
force of a classical Ohrnstein-Uhlenbeck type process. We 
remark that the first diffusion term turns out to be neg- 
ligible in many cases. For a discussion of that latter as 
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W(q,p) 




-0.5 



FIG. 2: Wigner function W(q,p) of a superposition of two 
Gaussian wave packets as in Fig. 1, here with qi — —5 = —§2, 
Pi = P2 = 0. The mutual distance is larger than individual 
spread. The oscillating pattern in between the two Gaussians 
indicates their coherence. 



well as various other limits of physical relevance we refer 
the reader to 10 . 



We may solve equation (3.12) with the help of the char- 
acteristic function or Fourier transform of the Wigner 
function, 



i{liq-vp)/h 



(3.13) 



Interestingly, diffusive effects on W can be confined to a 
Gaussian factor in x, 



Xt(fJ-,v) = (t(p,v)exp{ b(s)ds) (3.14) 
Jo 

x exp { a t^ 2 + 2f3 t fiv + 7t^ 2 }^ ; 

the cofactor ( then obeys the "Liouville" equation 
Ct = (—?- {ji/M + M + ^(Mil 2 + at)A (t , 



dv dfi 

(3.15) 

provided t he ti me dependent coefficients a, 0, 7 in the 
Gaussian (3.14) satisfy the linear equations 







{ 


-I 








-2/M 
-h 



2(MfiHo f ) -26 




For the transformation to ensure Coif 1 ^ v ) — Xo(a*j z/ )> the 
initial condition for these coefficients is (a , /3oj7o) = 
(0,0,0) at t = 0. 



We may solve the "Liouville equation" ( |3 .15 ) through 
the characteristic equations 



-(Am 2 + a t ) 
l/M b t 



whose integral yields the linear mapping 



M, 



(3.17) 



(3.18) 



clearly, the 2x2 matrix M* originates from the 2x2 
identity, M = ( 01 )» at t = 0. One can further estab- 
lish the identity det M t — exp ^ b(s)ds*j which will be 
useful below. 

Given an arbitrary initial density operator and thus 
the corresponding initial characteristic function xo(a*j v)-, 
we obtain the evolved x*(/x, v) from (3.14), using 



dfi a I dvoxo(jMi,i^o)8(jj, — fH)S(v — v t ), 

(3.19) 



with (fj, t ,i> t ) the solutions (3.18) of the characteristic 
equations ( [3.17 ). 



C. Coherence Norm for Distinct Wave Wackets 

The Wigner function W as well as its Fourier transform 
X c an b e employed to express the norms iVy introduced 
in (U) and [SHB] as 

N tJ {t) = (2nfi) [ dqdp\W^{q,p)\ 2 (3.20) 



= (2nh) / dfidvlxl'i^v)] 2 . 

The first of these expressions nicely shows that N\2(t) 
is indeed a good indicator for the appearance of coher- 
ences between wave packets, as it measures the weight 
of the absolute square of the oscillating pattern of the 
Wigner function in between the wave packets as shown 
in Fig. 2. The second expression is most convenient to 
actually ev aluate iVi2(t) for the oscillator model . In fact, 
the ansatz (3.14) and the general solution ( 3.19] ) yield 



N u (t) = (27r?i)exp{/ b(s)ds} (3.21) 
Jo 



with (fi t , vt) the trajectories from ( 3.1 8| ) 



With P l y 2 s (0) = I 



and representing the 
Gaussian wave packets (2.2), simple Gaussian integrals 
give 

IXo 2 (^o,^o)| 2 = (^TTh)- 2 x (3.22) 
exp{-(^ - d Q ) 2 /Aa - <r(ji - d P ) 2 /h 2 } , 
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where dq — \q% — q\\ and dp — \p% — Pi\ denote the 
separations in position and momentum of the two wave 
packets. We thus see the characteristic function of the 
coherence to be strongly peaked near those distances, 
with both widths of order ^/h (Recall that we had chosen 
minimum uncertainty wave packets with a ~ h). 



Evaluating the Gaussian integral in (3.21) results in 
the appealing form 



v l2 i/) = P(i)v^\-^{dp,d Q )Aj' 1 '' 



(3.23) 



for the coherence norm, revealing the quadratic depen- 
dence of the decay on the initial separations dQ,dp. The 
time dependence of this quadratic form is captured in the 
matrix At, which turns out bulky in this general case; it 
reads 



A t = C t 



with the matrix 



1/cr 
Aa/h 2 



Ct 



(3.24) 



C, = M? 



a t fit 
fit It 



(3.25) 



and M t the propagat ing m atrix from ( |3.1g| ). The func- 
tions (at,/3t,jt) hi (3.25) are the solutions of ( 3.16| ) 
with vanishing initial values (recall that the evolution 
of (at, fit, "ft) involves the diffusion coefficients Ct and dt 
while only the drift coefficients at and bt enter the deter- 
ministic equation for M t ). 

A prefactor 



V(t) = det 1 



1/cr 
Aa/h 2 



exp{ / b(s)ds} 
Jo 

(3.26) 

appears in ( |3.23| ) which is independent of the separations 
and to be identified with the similar expression (3.14- 
SHB) in [SHB], resulting from the Gaussian integration. 
As explained in detail in [SHB] , as long as the separation 
of the two wave packets is large compared to their indi- 
vidual spread, this V(t) may safely be replaced by unity 
for those short times for which N12 decays to esse ntially 
zero due to the relevant exponential term in (3.23). 

With the dependence of the coherence norm on the 
initial position and momentum separations dQ,dp now 
made explicit, we proceed to unveil the time dependence 
of the matrix A t in the exponent of the coherence norm 
( |3.23 ), both for the interaction dominated early-time 
limit as well as for the long-time golden-rule limit. We 
stress again that due to the large sepa rations (dQ,dp) 
we focus on, the prefactor V(t) in (3.23) can and will be 
replaced by unity in what follows. 



IV. LIMITING CASES 



A. Golden Rule: r sys < r dcc < r diss 



Our exact result ( 3.23 ) for the coherence norm allows 
us to investigate the decay of coherences on all time 
scales. Matters simplify considerably in the case of weak 
coupling, where decoherence between the two wave pack- 
ets may again be evaluated analytically. 

Let us first determine the coefficients (at, fit,^ft)- To 
that end we first recall that these coefficients vanish ini- 
tially and then realize that the inhomogeneities Ct,dt are 
of second order in the interaction. We may therefore re- 
place the 3 x 3 propagator matrix in the evolution equa- 
tion (3.16) in zeroth order, i.e. by entirely neglecting 
at and bt therein. The propagator matrix is then easily 
exponentiated to give 



ott 



fit 



It 



ds 



(c(s) (sin 2fi(t- s))/MQ (4.1) 
+d(s)(l - cos2fl{t - s))/M 2 fi 2 ), 
Is (c(s) cos 20(t — s) 



+d(s)(sin2fi(i - s))/MQ 
ds (c(s)MCl sin 2Q(t - s) 
-d(s)(l + cos2Cl(t - s)) 



Moreover, we may replace the t ime dependent coefficients 
of the exact master equation (3.9) by their lowest-order 
expressions, which are 



ds U[B(s),B]) cos n s +0(B 4 ) (4.2) 



Mfl 



1 



d t = 



ds(jr[B(s),B]) sin fls +0(B 4 ) 
ds({B(s),B})8ians +C(B 4 ) 
0(B 4 ). 



2MVL 
t 

ds±({B{s) 7 B}} cos Os 



Note that perturbati on t heory preserves the correct 
short-time expansion (3.1C). 

Now t hat the coefficients at, fit, It within the matrix Ct 
in ( 3.25| ) are revealed as of second order in the interaction, 
it suffices to replace M t by its zeroth-order approximant 
in that exponent. Neglecting the second-order terms at 
and bt in the definition ( 3.1 8] ) , we find 



Mi 



cosfH — AffisinfH 
sinfM/Mf2 cosfii 



(4.3) 



With (4.1) and (|4.3|) we have access to the full time de- 



pendence of the decohering quadratic form in the expo- 
nent of the coherence norm ( 3.23| ). 
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Of particular interest is the long-time limit, ilt ^> 1, 
which allows for many system oscillations up to the ob- 
servation time. In that limit, furthermore assuming the 
bath correlation time to be much shorter than the system 
time scale fl~ , we may safely replace the time depen- 
dent coefficients c(s) and d(s) by t heir asymptotic values 
c(oo), d(oo) under the integrals in (4.1). Moreover, we see 
from the latter expressions that for flf ^> 1 all oscillating 
terms vanish and the only surviving contribution arises 
from the constant terms involving d(oo). The Golden- 
Rule limit fit ^> 1 thus finally yields the coefficients 



at 



d(oo) 

M 2 n 2 
0t = o, 
7t = d(ao)t 



t, 



(4.4) 



and the matrix ( |3.25 ) 

C t = d(oo)t 



1 



(4.5) 



With Ct alre ady o f second order in the in the interaction, 
we see from ( 3.24 ) that A t coincides with Ct in this or- 
der of perturbation theory. The well known exponential 
decay 



Ni 2 (t) = exp 



d(oo)t 



d 2 Q +d 2 P /M 2 ft : 



-t/r 



(4.6) 

characteristic of golden-rule decoherence results (recall 
that we may drop the slow ly varying prefactor V(t) of 
the general expression (p3|)). Clearly, since we allow the 
system to evolve for many cycles, any distinction of the 
position as the coupling agent in the interaction Hamil- 
tonian has disappeared. Through the sequence of system 
cycles, position and momentum interchange their role pe- 
riodically such that in the long-time limit of decoherence 
only the joint quantity ((c£q) 2 + (dp) 2 / (M 2 Vt 2 )) appears 
as acceleration factor. 

For the final result for the golden-rule decoherence time 
we still have to determine the rate d(oo). In the relevant 
order of perturbatio n theor y and upon using the bath 
correlation function ([3.4[|3.5|) we get 



d(oo) 



ds U{B(s),B}} cos tts 



= hM 1 nf c (n/A){n th + ^): 



(4.7) 



where again nth = (e nn ^ kT — l) -1 is the thermal number 
of qu anta in the oscillator and f c (il/A) — » 1 according to 
(3.6) . Comparing with the lowest-order drift coefficient 



6(oo) - 



Mfl 



ds-([B(s),B}) sin Ms 



= 7 / c (fi/A) -> 7 



(4.8) 



we recover the previously announced interpretation of 7 
as the classical damping constant, 7 = 1/Tdi ss , as well as 



the well known golden-rule decoherence time already 
presented in (1.2-SHB) of [SHB]. That latter expression 
implies the familiar golden-rule time-scale ratio for deco- 
herence and dissipation, 



'dec 
GR 
diss 



des 



(4.9) 



we see decoherence accelerated over dissipation by 
th e squared ratio o f a de Broglie wavelength A t h = 
y / h/MH,(n, t h + 1/2) and an effective distance between 

the two wave packets, c? G ff = ^Jdg + d 2 p /AI 2 Q 2 . 



B. Interaction dominance 1: 



Tdec Tsys , 7"dis 



As soon as the initial separation (dq, dp) between wave 
packets extends beyond quantum scales, decoherence is 
dominated by the interaction Hamiltonian Hi nt since the 
free motion generated by -ff S ys eventually becomes negli- 
gibly slow by comparison. This is the regime we discussed 
at length and in a general setting in the accompanying 
paper [SHB]. We must recover those previous findings 
for the exactly solvable oscillator mod el by suitably sim- 
plifying the general expression (3.23) for the coherence 



During these very short initial time spans, we take into 
account dynamics on cl assical time scales due to the "Li- 
ouville equation" ( 3.1 5| ) to lowest order in t only a nd rc - 
place the corresponding propagating matrix M t in (3. IS) 

by M t = ( t j M ~ M ^ 2t \ + Q{t 2 ). Next, we have to deter- 
mine the coefficients at,/3t and 7 t in ( |3.16 ). As before, 
we expand all dynamical quantities connected to the os- 
cillator dynamics to the lowest relevant order in t and 
find 



a 1 



0t 



It 



ds a{B(s),B(0)})(2t + «)(*- s) 2 /3M 2 , 



ds a{B(s),B(0)})t(t~s)/M, 



- / ds({B(s),B(0)}). 
Jo 



(4.10) 



We see the importance of the details of the (real part of 
the) bath correlation function, (^{B(s), B(0)}), describ- 
ing dynamics on reservoir time scales. 

The matrix ( 3.25| ) determining the decay of coherences 

is 



ds a{B(s),B(0)})(t-s) 



(t- s)(2t + s)/3M t/M 



(4.11) 



t/M 



As in the weak coupling case, the difference between 
the relevant matrix A t in (3.24) and Ct is negligible 
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for the fast decoherence resulting fro m lar ge separations 
{dp, d p), and with A t replaced by ( 4.11 ), we get from 
( p3| ) the cohe rence norm as the exponential 



N 12 {t) 



exp 



d Q d P t/M 



ds({B(s),B(0)})(t-s) (4.12) 



d P 2 (t 



s)(2t + s)/6Af 2 )}, 



see [jTlJ and [SHB] for the case dp — 0. As before, we 
may neglect the s lowly varying prefactor V{t) of the exact 
expression ( |3.23 ). 

Two remarks about this expression are in order. First, 
the oscillator model clearly confirms the general result 
(6.5-SHB) (for dp = 0) of the accompanying paper 
[SHB] , that derivation being based solely on the assump- 
tion Tdoc Tsys;7"diss and the Gaussian character of the 



bath coupling agent. Secondly, ( 4.12| ) may be recog- 
nized to coincide with the squared absolute norm of the 
Feynman- Vernon influence functional |?], with the 
classical paths represented by the short-time expression 
qt = qi + p\t/M (replace one by two for th e seco nd path 
q' t ). Examples of decoherence following ( 4.12 ) will be 
shown in Sect. VI. 



Interaction dominance 2: 



, Tdis 



In the extreme case when decoherence is even faster 
than any environmental time scale it is sufficie nt to ex- 
pand the whole matrix A t in the exponent of (3.23) in 
powers of the elapsed time t, mirroring the corresponding 
expansion of the logarithm of the interaction propagator 
in [SHB]. Crucially, we replace the correlation function 
i({B(s), B(0)}) under the integral by its initial value 
(B 2 ). Starting from (4.11) of the last section we find 



C t = (B 




(4.13) 



in each entry neglecting higher order terms. Again, the 
short time approximation demands A t to be identical to 
B t from ([4.13) to the relevant order, and the general 



result (3.22) turns into 



N 12 (t) = exp|-^(^) (4.14) 
(d Q 2 t 2 + d Q d P t 3 /M + d P 2 t 4 /4:M 2 ) } 



which is the universal law (3.14-SHB) for the decay of co- 
herences we found under very general conditions in the 
accompanying paper [SHB] and for dp = in jtl). We 
refrain from including the prefactor V(t) here, that fac- 
tor in (3.14-SHB) being safely replaced by unity for the 
relevant short times characteristic of the decay for asymp- 
totically classical initial separations c?q, dp. 

Any coherences present in the system will have van- 
ished according to (4.14), before the quantum state can 



be aware of any potential V(Q) - this is why the oscilla- 
tor frequency f2 has disappeared entirely from the final 
expression (4.14) (and also from ( 4.12| )). Recall that sys- 
tem oscillations with fre quenc y Q were crucial for the 
golden-rule result. Eq. ( |4.14 ) also confirms the differ- 



ent scalings of the corresponding decoherence times r ( 



Q 

dec ' 



'dec*' an d r dcc discussed in great detail in the fourth sec- 



tion of [SHB]. 



V. AN ALTERNATIVE MEASURE FOR 
DECOHERENCE 

Before turning to actual examples, we replace the gen- 
eral coherence norm N12 (t) by a somewhat simpler quan- 
tity with equal power to reveal decoherence. 

As we are going to restrict ourselves to superpositions 
of symmetrically located wave packets with q 2 = —qi, 
and p 2 — —pi, an alternative indicator of the decay of 
coherences is the value of the Wigner function of p sys at 
the origin, 



Wt(0,0) 



1 

itTi 



dq(q\p\-q) = 



1 



d P {p\p\-p). (5.1) 



Apparently, Wt(0,0) is a measure for the weight off the 
diagonal, both in the position and momentum represen- 
tations; it is linear in the density operator and hence 
easier to access numerically, given the master equation; 
it is intimately related to the previously employed N12 (t) 
as we will briefly show. By definition, we have Wt(0, 0) = 

/ dixdvxtiw) = 2k I d ^ diy J2i 3 (HCjXt (M,^)- Thc 
diagonal terms x} 1 an d Xt 2 do not noticeably con- 
tribute to W t (0,0) since for superpositions of far-apart 
wave packets their Wigner correspondant is located near 
(qi, pi), and (— qi,— pi), respectively, i.e. far away from 
the phase space origin, as is also apparent in Fig. 2. 
Furthermore, for the symmetric case considered here, we 
find real Gaussians for the characteristic functions such 
that Xt 2 = X? 1 = \x 21 \ an d therefore 



Wi(0,0)oc / dpdv\ X \ 2 {^v)\. 



(5.2) 



Comparing the latter with expression ( 3.20 ) for the 
norms Nij we see that for the symmetric superposition of 
wave packets considered here, the linear quantity VFt(0, 0) 
has equal power to reveal the fate of coherence as the co- 
herence norm Ni 2 (t). Going through the very same steps 
as in the previous section we find 

(5.3) 



m 2 (t) = W t (o,o)/Wo(o,o) 
1 



V{t) exp 



2 {d P ,d Q )A t 



dp 
d Q 



with the sa me p refactor V{t) and matrix At as in the 
expression (3.23) for Ni 2 (t). Neglecting the irrelevant 
prefactor V(t) reveals that the simpler quantity n\2{t) 
is essentially just the square root of the coherence norm 
Nx^if) and our previous discussion of that quantity in 
Sect. Ill equally applies to ni2(i). 
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VI. NUMERICAL RESULTS FOR THE 
OSCILLATOR MODEL 



We illustrate numerically the decoherence of a super- 
position of two symmetrically located coherent states, 
= | a) and <f2 = \ — a), which are labeled by a com- 
plex number a, | a) = e~2l Q l e aa |0), see M. The dimen- 
sionless creation and annihilation operators, for instance 

a = (^/MQ/2h Q + iP/V2MQhj , contain the quantum 

reference scales of length and momentum characteristic 
of the harmonic oscillator. We may employ those scales 
to introduce dimensionless distances between the two co- 
herent states 



d Q = dq/y/h/Am = V2{a + a*) 
dp = dp/VMMl = V2(a - a*)/i , 



(6.1) 



which will be the relevant quantities in the following. In 
these units, large separations with respect to the spread 
of coh erent states, as required for our general formula 
( 3.23 ), simply means large compared to unity. Truly 
macroscopic separations could lead to dimensionless sep- 
arations of the order, say, dQ ~ 10 17 , a scale way beyond 
our examples below. We shall let dQ or dp range from 
16 to 4000 and in one extreme case to 800000. It will 
become clear that the transition from slow golden-rule 
decoherence to rapid interaction dominated decoherence 
is fully captured in the mesoscopic range considered. 

We assume a thermal environment with all parame- 
ters of the total Hamiltonian fixed once and for all. It 
is the initial state of the system oscillator only that we 
alter between simulations. The transition we are after 
is traced out by varying the dimensionless initial phase 
space distances (dQ,dp) between the superposed wave 
packets. 

All simulations are performed with the small damping 
rate 7 = 10 _5 17, a rather large temperature according 
to ksT = 20hn, and an environmental cutoff frequency 
A = 100f2. The latter is the largest frequency involved 
in either system or environment. We expect interaction 
dominated decoherence according to (4.12) to set in as 
soon as the decoherence time scale becomes shorter than 
the shortest system time scale, here O -1 . The extreme 



case of quadratic decay according to (4.14) (or qubic or 
quartic, depending on the nature of the superposition - 
see also [SHB] ) occurs as soon as Td ec is even shorter than 
reservoir time scales, here determined by A . The nu- 
merical solution of the exact master equation is based 
on the corresponding exact stochastic Schrodinger equa- 
tion for this model, as established recently |Oj. This ap- 
proach has the nice feature of preserving Gaussian pure 
states within each run of the Schrodinger equation for 
a given realization of the driving noise and, moreover, 
allows easy access to the otherwise involved expressions 
for the time dependent coefficients entering the evolution 
equations. 



A. Golden rule 

We can estimate for which dimensionless effective sep- 
aration d c fi = \J 3q + dp the golden rule result should 
apply for our choice of parameters. ;.From (|4 

(O) we find = ( 7 (n th + l/2)<4 

the golden rule assumes a long-time limit fit ^> 1. We 
can thus trust the golden rule decay law only as long as 
^ T dcx 1 which puts an upper bound on d e g. Together 
with the lower bound arising from the requirement of a 
large separation we have 



1 um ^.uj ) and 

, \ -1 

2 ' Recall that 



1 < d cff < v^TKh + 1/2). 



(6.2) 



For the choice of parameters above, we find 1 <C d e g -C 
70, a fairly narrow band for the applicability of the golden 
rule. 
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FIG. 3: Golden-rule decay of an initial superposition of 
two coherent states. Distances are dQ = 16 and dp — 0. 
Environmental parameters are 7 = 10 _5 17, A = 10017, and 
kT — 20ftl7. Full line ("ex"): exact numerical calc ulation, 
indistinguishable from the analytical expression (3.23) (dash- 
dotted, "an"). Dotted line ("GR"): golden-rule result (4.6) 
which fits smoothly through the slightly oscillatory exact 
result. Interaction-dominance results (dashed lines, "idl", 
"id2" ) are irrelevant here; the influence of the system Hamil- 
tonian is crucial. 



In Fig. 3 we show the decay of an initial superposition 
of two coherent states with separations rf e ff = dQ = 16 
and thus dp = 0. The full line ("ex") shows the ex- 
act numerical result. Apart from some tiny oscillations 
around the exponential decay the golden-ru le ex pectation 
(dotted and labeled "GR") according to (4.6) is nicely 
confirmed. Note how the tiny oscillations loose final rel- 
evance as soon as fit ^> 1. These oscillations have their 
origin in the periodic change of the system state between 
a superposition with respect to Q to one with respect to 
P. Decoherence is more efficient for Q-superpositions, 
and less effective in those time spans when the system 
state assumes a superposition with respect to P. Dur- 
ing those periods decoherence is slowed down slightly. 
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During one oscillator period there are two swaps be- 
tween a Q- and a P-superposition, which fact explains 
the frequency 2fl of the oscillations in question. The 
dash- dotted line ( "an" ) depicts the analytical expression 
(3.23) and is almost indistinguishable from the exact re- 
sult. The two dashed curves ("idl", "id2") correspond to 
predictions from the two interaction dominated regimes. 
They are clearly irrelevant for the still "slow" decoher- 
ence of the underlying "microscopic" superposition. The 
system Hamiltonian has a strong influence on the decay 
in that regime. 



efficient again after another quarter of the period, and 
periodically thereafter. The periodic alteration between 
slowdown and acceleration is now a prominent feature 
rather than a small fluctuation. For times short com- 
pared to the system time scale, fit <C 1, the influence of 
the system Hamiltonian is negligible and thus the gen- 
eral result (4.12) of interaction dominance (dashed line, 
"idl") initially applies as can be seen from the good 
agreement for short times. 



Interaction dominated decoherence 1 



B. Crossover regime: Tdcc ~ T sys 

The golden rule ceases to be relevant as soon as the 
initial separation between the superposed wave packets 
increases to such an extent that the decoherence time 
scale Tdcc becomes comparable to r sys . Then the decay 
must still be system specific, that is depend on the de- 
tails of the dynamics induced by the system Hamiltonian 
and the time scales emerging from the bath correlation 
function. It is the special virtue of the oscillator model 
to allow rigorous treatment of this crossover case. 




FIG. 4: Crossover regime, dg = 80 and dp = 0, otherwise 
same parameters as in Fig. 3. Golden rule no longer valid, uni- 
versal interaction dominated regime not yet reached. Full line 
("ex")for exact numeric al ca lculation; indistinguishable dash- 
dotted line ("an") fo r (|3.23| ). Only for small times, fit <^ 1, 

becomes visible. 



decay according to (4.12) ("idl" 
tions due to system Hamiltonian, see main text. 



Oscilla- 



Fig. 4 shows the decay of an initial superposition of 
coherent states with an intermediate position separation 
cLq = 80 and equal momenta, dp = 0. The (indistin- 
guishable) full and dash-dotted lines show the exact mas- 
ter equation result and the evaluation of the exact ana- 
lytical expression ( 3.23 ) . We see that decoherence due to 
coupling to the position is very effective initially; how- 
ever, as the system dynamics moves the Q-superposition 
to one with respect to P after a quarter of the oscilla- 
tor period Sit w 7r/2, decoherence slows down. Turning 
the state further in phase space, decoherence becomes 



As soon as the separations d,Q or dp between the su- 
perposed coherent states are sufficiently large, the system 
Hamiltonian becomes irrelevant and interaction domi- 
nates, see Fig. 5. For a separation of dq = 200 and 
equal momenta dp = 0, decoherence is faster than any 
system time scale, yet still longer than the environmental 
correlation time A -1 . This is the regime that can fully 
be described by the general decay l aw ( 6.5-[SHBl) of the 
accompanying paper [SHB], here ( 4.12 ). As mentioned 
in [SHB] , unless the spectral density differs fro m zer o at 
zero frequency - which is not the case here, ( 4.12 ) de- 
scribes non-exponential decay. We further remark that 
the golden rule (dotted line, "GR" ) now predicts far too 
slow decoherence. 




FIG. 5: Interaction dominance, d<g = 200 and dp — 0, 
otherwise same as in Fig. 3; note change of time scale. Full 
line ("ex") for exact numeri cal ca lculation, indistinguishable 
dash-dotted line ("an") for (3.23). Golden rule (dotted line, 



"GR") now meaningless; interaction dominance (dashed line, 
"idl") thro ughou t whole relevant time span. Extreme short- 
time result (4.14) (dashed line, "id2") enters validity for very 
early times. 



D. Interaction dominated decoherence 2 

We expect univer sal d ecay according to our simple an- 
alytical expression (4.14) (corresponding to result (3.14- 
[SHB]) in [SHB]) as soon as the decoherence time is the 
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shortest time scale involved, here even shorter than the 
inverse environmental cutoff frequency A -1 . We thus find 
lower bounds for the size of superpositions that indicate 
the beginning of this universal interaction dominated de- 
coherence regim e. 

We use (3_/7) to find an expression for (B 2 ) and con- 
clude from the discussion of interaction dominated deco- 
herence times in [SHB] that for a superposition in po- 
sition space, Ar^ ec — \JhMl/ (kTj) /d~Q . The latter has 
to be small compared to one in order for this universal 
regime to be important. With the choice of parameters as 
before, we find the condition d,Q ^S> 700 for the universal 
Gaussian decay law to be relevant. 
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0.004 0.006 

at 



0.008 



0.01 



FIG. 6: Interaction dominated universal decoherence, <1q — 
4000 and dp = 0, otherwise as in Fig. 3; note change of time 
scale. Four results in agreement: full line ("ex") for exact 
numerica l calc ulation, dash-dotted line ("an") for analytical 
formula (3.23), two dashed lines ("idl" and "id2") for inter- 
action dominance, the latter the Gaussian exp{— (i/ r <£ c ) 2 } 
from (4.14). Golden rule (dotted line, "GR") grossly under- 
estimates life time of coherence. 



In Fig. 6 we show the decay of an initial superposition 
of coherent states with a yet larger separation cLq = 4000 
and dp — 0. Note the difference in time scales along 
the i-axis between the various figures. While the full 
line ("ex") represents the exact numerical calcul ation , 
the dash-dotted line ("an") is our exact formula ( p. 23 ). 
In this case, interaction dominance holds (dashed lines, 
"idl" and "id2"), the extreme short-time result Q4.14 ) 
suffices to describe the decay over the full relevant time 
interval. The latter is a Gaussian exp{— (£/ T dL) 2 } decay 
law. Clearly, the golden rule is not applicable here and 
wrongly predicts far too fast decoherence. This is due to 
the quadratic golden rule scaling of the decoherence times 
with distance, as compared to only a linear dependence of 
the true interaction dominated decoherence time [SHB] . 

A superposition in momentum space with equal posi- 
tions d,Q = is more stable under the position dependent 
coupling H mt = QB between system and environment. 
To reveal interaction dominance and our universal decay 
law, a momentum superposition will have to be much 
larger (in our dimensionless units) than a position space 



superposition. As in the previous case, we get an esti- 
mate for the border of applicability of our universal law 
(see the discussion in the accompanying paper [SHB]) 

At[ cc ~ (M 3 /(fc™ 7 )) 3 /y/fo. This quantity has to be 
small compared to one for universal "quartic" decay. For 
our choice of parameters, dp >• 70000 is required. 



-nO.5 




0.002 



0.004 0.006 

£lt 



0.008 



0.01 



FIG. 7: Interaction dominated universal decay of initial su- 
perposition of two coherent states with opposite momenta and 
equal positions. Distances are dQ = and dp = 800000, other 
parameters as in Fig. 3; note change of time scale. Four re- 
sults in agreement: full line ("ex") for exact numerical calcu- 
lation, dash-dotted line ("an") for analytical formula ( 3.23 ), 
and a dashed line ("id2") of interaction dominance, t he lat - 



ter now the quartic exponential exp{ — (t/ T dD } l aw (1-14) 



Golden rule (dotted line "GR" ) underestimates lifetime of co- 
herence even more drastically as in Fig. 6. 

We check on that case in Fig. 7 which depicts the de- 
cay of an initial superposition of two coherent states with 
momentum separation dp = 800000 and equal positions 
dQ = 0. As before, the full l ine (" ex") is the exact numer- 
ical result compared with ( 3.23 ) (dash-dotted line) and 
the interaction dominated results (dashed line, "id2"). 
Again, the golden rule by far underestimates the deco- 
herence time, due to the wrong scaling: while scales 
with dp, the true T^ ec , as it may be determined from our 
approach to interaction dominated decoherence [SHB], 

scales with the square root \J dp only. However, Fig. 7 
does confirm the quartic deca y n 12 (i) = exp{— (i/ T doc) 4 } 



we predict according to ( 4.14 ) 



VII. CONCLUSIONS 

The exactly soluble model of a damped harmonic os- 
cillator coupled to bath of oscillators illustrates our pre- 
diction of various decoherence scenarios. For microscopic 
distances between the superposed wave packets, decoher- 
ence is so slow as to witness many free-system cycles and 
follows the standard exponential golden-rule decay. As 
the separation increases, decoherence and system time 
scales eventually become comparable and a more compli- 
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cated, system-specific decay of coherences results. Our 
findings indicate that golden-rule predictions may fail for 
even only moderately sized superpositions. For further 
increased separations we enter interaction dominated de- 
coherence as presented in |n| and in much more detail 
in the accompanying paper [SHB]. Now Td cc *C T sys and 
decoherenc e pro ceeds according to a system-independent 
decay law (4.12), incorporating the bath correlation func- 
tion. In the extreme case when decoherence even outruns 
intra-bath processes, the decay becomes a simple Gaus- 
sian for an initial separation with respect to position, and 
an initially slower fourth-order exponential if momentum 
is the differing property; position is distinguished over 
momentum since it is taken as the system coupling agent 
in the interaction with the bath. 

The special virtue of the oscillator model studied here 
is to allow inspection of the crossover between the long- 
time golden rule regime valid for microscopic superposi- 
tions, and the interaction dominated regime relevant for 
more and more macroscopic superpositions. 

As experimental efforts towards resolving decoherence 



dynamics continue H], it will be fascinating to see 
whether larger-scale superpositions can be realized and 
whether decoherence dynamics can be pushed towards 
the interaction dominated regime illustrated here and in 
[SHB]. 
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